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Abstract 

Bose gas in a random external field is considered. The sigma model like effective action both for 
weak and strong random fields compared with the interaction between particles is derived by aver- 
aging over the random field and integration over the quantum fluctuations in the framework of the 
Keldysh-Schwinger technique for irreversible processes. Using this effective action the properties 
of the system in strongly fluctuating random field are analyzed. 







I. INTRODUCTION 



The problem of transition between localized and extended states of many particle systems 
attracts attention of condensed matter theorists and experimentalists for a long time. The 
theory of non interacting electrons in the random external field in one and two dimensions 
p.edicts localization otpa^cte fa a.Mt.„ny weak d.o.de. g g g In th.ee dimension, 
the state of non interacting electrons is localized for sufficiently strong disorder [5|], and 
the state is extended for the disorder less than a critical edge. For negative energy with 
absolute value much larger than the mobility edge, the system is strongly localized and 



its density of states is described by so called "tail" of density of states 



fl, S, B, fl. 



The common point of view is that the repulsive interaction between electrons suppresses 
ocalization. So far this problem was mainly considered in the case of week disorder [12I . 
I3I . However the problem of investigation of transition between localized and extended 
states caused by interaction between particles in three dimension is of great importance and 
interest, but it cannot be solved in the framework of the weak disorder. 

The problem of localization in a random field was investigated for the case of the Fermi 
statistics of particles. In this case if the Fermi energy is large enough the interaction between 
particles can be neglected, as it was done in[l|, I2I, js], [ol, 10], 11]. The analysis of the 
problem when both random field and the interaction between particles are involved, is the 
sufficiently more complicated. 

In this paper we consider the case of Bose statistics of particles, although, the technique 
developed here, in general points, can be applied for the Fermi case too. For Bose system 
in random field the interaction between particles plays fundamental role and cannot be 
neglected. This is because the nonzero energy of the Bose system without random potential 
is defined by the interaction between particles, and exactly this interaction prevents the 
fall of the system to the very bottom of the fiuctuating random potential wells. Our goal 
is development of a technique for the description of the transition from the localized state 
to the delocalized superfiuid state. The experimental and theoretical investigations of this 
Tansition in the systems of different nature are of an interest for a long time 28| - 36| . 



|2l|, 
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2J], [25|, |23|, |26|, |27|]. In the works [22|] - [23| the Bose gas in a small random 
external field was considered and corrections to the temperature of the superfiuid transition 
and to the superfiuid density were calculated using perturbation theory in random field. In 
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the recent works 26|, 27|] the case of the strong random field, or small density of Bose gas, 
were considered quahtatively. The technique developed here enables to consider both strong 
and weak random fields compared to the strength of interaction between the Bose particles. 
The idea to average over the random field at the first step of calculations if the Keldysh- 
Schwinger technique [3], jl5| was already exploited previously We develop this 

idea and obtain the effective action in the form of a sigma model, which describes both the 
cases of the weak and strong random field. In this work we mainly focus on the application 
of the obtained effective action to the case of the strong random field. 



II. AVERAGED GENERATING FUNCTIONAL 



The Keldysh-Schwinger technique 151], [l6|], along with supersymmetry and replica tech- 
nique, is known as an effective in the problems of localization in a random static potential 



171], 18|. In this technique absence of denominator in the representation of generating func- 



tional makes it easy to take an average over the Gaussian random static external potential 
V {~r^), provided some natural assumptions are made, like absence of interactions for asymp- 
totic in-states which are turned on adiabatically. The generating functional for the system 
in random static potential V (T^) is 



Z [J, V] = j D0D0exp {iS [V] + iSj} 



with the action 



S [V] = So + ASo [V] + S,nt 
and interaction with the sources Sj. Here 5*0 is the free particle action 



5n 



dtd'^r 



ASq [V] is the interaction of the particles with a random external field 



ASo [V] 



dtd'^r {(pV(j)} 



(1) 
(2) 

(3) 
(4) 



and Sint and Sj are the interaction between the particles and their interaction with the 
external source fields, which are auxiliary and put to zero at the end of calculations. 



Sint = - \gjdtd''r {(t){t,-r)(l){t,l^)(l) (t, r ) (t, r ) } 



(5) 



Sj = jdt(fr [J0 + J] (6) 

In ([3])- ([6]) the time t integration contour is the standard Keldysh-Schwinger double time 
contour with the return, going from t = —oo to t = +oo (shifted up) and back to t = —oo 
(shifted down) in the complex t plane. The complex particle fields (p, (j) and their sources 
J, J are functions of the time and space coordinates t, ~r . At the moment we don't fix the 
statistics of the particles, but in what follows we concentrate on the Bose particles case, g 
is coupling constant, and the function ^ ^"p^j is the excitation spectrum of the ideal gas, 
^ ("^) ~ 2rn ~ where the momentum operator is "p^ = — zV , and /i is the chemical 
potential. We choose the system of units in which the Plank constant h and the mass of 
particles m are equal to unity h = m = 1. 

The probability of the static random field distribution is supposed to be Gaussian with 
the correlator of white-noise type 

<V{1^)V (7") >= (l/Ko) - 7") 

with some constant Kq: 



P [V] = ^ S T (7) 



exp 1 






\ 


/ DVexp 




1 d'^rd'^r' {V {'r)f^ 



As we said, the absence of denominator in ([T]) makes it possible to integrate over the random 
potential V in the averaged generating functional: 

< Z [J, V] >v= j D^D'^DVexp {iS [V] + iSj} P[V]= j D(j)D'^exp {iS} (8) 

where 



S — So + Sint + Sav + Sj 



(9) 



The parts of the action 5*0, Smt, Sj are given by equations ([3]), ([5]), ([6]), and new additional 
interaction term Sav generated by random field averaging, has the form 



Sa 



2Ko 



d\ 



<j)dt(j) (t, ~r) (j) (t, ~r 



dt'(f) {t\ V) (f) {t\ ~r 



(10) 



It is convenient to rewrite the action 5* in the so called three-angle representation 16| . 
[lol . and introduce the fields \1/ and ip instead of the fields 0i and 02 , where 1, 2 are the 
indexes of the upper and lower brunches of the time contour respectively: 



J 



+ 



2 ' 
Ji + J2 
2 ' 



^ = 01 - 02 
j = jl - 32 



In these variables the parts of the action take the form 



Si 



Sr.. 



int 



2Kn 



d\ 



dt + < / dt' [-^ij + ip-^] 



Sj = j dtd'^r { {J^P + ^J) + (j^ + ^j) } 



(11) 



(12) 



(13) 



(14) 



(15) 



where the inverse retarded and advanced free field Green functions Gqh, Gqa in frequency 
representation are 



(poA 



UJ + fi + —V^ + n 
2m 

1 ^2 
to + fi+ - — V - 27 
2m 



(16) 



with infinitesimal positive 7. In time representation u should be replaced by idf. In the 
three-angle representation asymptotic of the fields ■0, ■0 at t — 00 satisfy the condition 



ip {t 00, 1^) = ip {t ^ 00, ~r) = 0, and the fields \1/ are non zero at t ^ cxo [17|. Note 
that the classical fields, i.e. the fields which in functional integral technique correspond not 
to operators but to c- number values, have ip = (pi — (p2 equal to zero at any time moment. 
The fields -0, ip describe quantum fiuctuations of the quantum system, and \1/ can be 
considered as the classical components of the fields 



18l | . In the case of large occupation 



of some degrees of freedom the quantum fiuctuations 0, %p corresponding to these degrees of 



freedom give small corrections compared to the classical components ^E', ^E', and in principle 
can be taken into account in the lower orders of perturbation theory. If the system has Bose 
condensate it is described by the fields \E'. 



III. INTEGRATION OVER QUANTUM FLUCTUATIONS. 

Application of the Keldysh-Schwinger technique for localization problem in Fermi systems 
in weak random field 18l was based on the construction of so called Q-Lagrangian l9| , like 
for example in replica or supersymmetry technique [9|]-^J]. Correctness of application of 
this technique to strong localization problem is doubtful. The main goal of this paper is to 
consider the properties of many-particle system in strong random field. Here we consider 
Bose gas of small density in strong random field. We suppose nq^« 1, where n is the 



averaged density of the gas, g is related to the scattering length [20|, and the fluctuations 
of random field are assumed to be larger or of the order of /i,. = ng. 

The novelty of our approach is that we start from integration over the quantum fields ip, 
ip, getting an effective action for the classical fields \E', \E'. In order to integrate over ip, ip we 
introduce new fields r], fj by 



(t, T^) = ^ (t, T") {t, V) (17) 

1 

i) (t, ~r) =r] {t, T") ^ (t, T") 



We preserved the form of equations which will be valid for both statistics of the fields ip,"^, 
although our goal in this paper is, as we said, the bosonic case. In this context the above 
equations are written in the way which avoids a division by a Grassmann variable in the 
fermionic case. Here and below the form of equation f|T71) is also convenient in the case of 
additional degrees of freedom for the field \E', like spin. The action Sav is Gaussian in t], 
r]. Our strategy will be to integrate over the fields 1], fj and write the effective action in 
terms of the chiral fields. Without the interaction term Smt this integration is exact both 
for Fermi and Bose cases. The effective action takes the form of a sigma model. In this 
paper the effective action will be obtained for the Bose gas with the repulsive interaction 
between particles, but the technique used here works in fermionic well. 
The Jacobian of transformations (fTTI) is 



In terms of new variables rj, rj the action Sav takes the form 



Sav = ^J d'^r J dtj dt' [r] {t, r)+r] {t, T")) {r) {f , + r] {f , T)) (19) 

This is the crucial point of the calculation: there are two factorized time integrations in the 
last equation. 

For future convenience we rewrite the action, introducing new (chiral) fields e, e, /, / 
instead of the classical fields \1/ : 

^ e = "*— ^ (20) 



/ = expQA); 7 = exp(^-iA) (21) 

A ^ In (-pj) (22) 

The constant independent of 7^, is introduced to make the argument of logarithm 
dimensionless. The chiral fields satisfy the following constraints 

ee = l, 7/ = l (23) 

At this stage we will neglect the highest order terms of the quantum fluctuations ip in Sint 
( fT3|l proportional to It is quite well understood [3], 1^, 13], [l8] that these terms are 
responsible for relaxation processes induced by interaction between the particles . These 
interaction responsible terms can be neglected if the density of the system or the interaction 
between particles are small ng^ << 1. Lets emphasis that this approximation, nevertheless, 
give the possibility to consider random field fluctuations of any strength exactly. They are 
essential if the kinetics of relaxation to an equilibrium state is considered. In this paper we 
suppose that the system is at an equilibrium state and there are no relaxation processes. 

With this assumption one can see that in terms of new fields (!20!) - (l22l) the free and the 
interaction parts of the action f|T2l) . f|T3l) combine to 



where 



So + S^nt = J dtd^r {fjL + Lt]} (24) 
L {t, r ) = 7 (tdj) + e (tdte) + + ( V ^e) + (25) 



L {t, y) = {tdtf) f + i-tdte) e + /i + ^ ( V'e) e+ (26) 



L + L = e (idte) + {-idte) e + 2/i - ^ Ve) Ve) + (27) 

+ i/V V - 2^7n(°) (/)^ 
m 

The constant n^c^ with dimension of density is arbitrary and actually drops out from the 
equation (1271) because of the gradient operator and Eq. (pT]) . We will fix its choice below. 

Eq. fl27|) was obtained from Eqs. fl25l [261) using the constraint (!23|) on the fields e, e. 
The remaining part of the action Sj can be written as 

Sj = J dtd'^r { {nfr^l^ [Jer] + rJeJ)7+ {n^c'^f'^f {je + ej) } (28) 
Equations (120|) can be written in polar coordinates p, (p for complex fields ^! 

= pexp (iip) = pe; \E' = pexp {—i^p) = pe; (29) 
p = ^/W^ 

Note that the functional integration measure changes under the transformation \l/ — ^ A, 

as 

Y[j^^D^D'^DT]Dr] = Y[D\D<pD7]Dr] 



If we split out the zero frequency mode (time independent) part rjo, Lq of rj, L 

Vit, r) = -Voir) + Sri{t, r); r]{t, r) = -r]^ {r) + 6r]{t, r) (30) 
L {t, r) = Lo (7") + 6L {t, r ) ; I (t, 7") = Iq (7^) + 6L {t, 
with the total time evolution T ^ oo, the action Sav (fT9l) becomes just 

Sa. = ^J d'r{r^,ir)f (31) 
whilst the other terms of the action take the form 

So + S,nt = j d'^r ||MC^±M j + y d^rdt {6t]5L + SLSr]) (32) 

Since the last term in ( |32l) is the only one where Srj, dfj appears, functional integration 
over them gives additional constraint SL = 5L = 0. It means that only time independent 
(zero frequency time Fourier components of the) functions L, L contribute to the functional 
integral. However, in general it doesn't mean that fields e, e, /, / are time independent. L, 
L can be taken at any time moment, e.g. at t = oo. With this choice we will omit the time 
argument for the functions e, e, /, / below. 

The first two terms of L + L Eq. ( 1271) in the first summand can be transformed to a renor- 
malization of the chemical potential: In terms of polar variables (l29l) they are transformed 
into 

i I d^dt [e{^^te) + {-id,e) e] = -2 | rf.,^(^^iZl_^iz^^^, i^) 

The difference of phases Ayj {~r) = ip (cxd, ~r) — Lp (— oo, 7*) can be written as Ayj (7^) = 
Sip (T') + 27rn, where < (5<y9 (T') < 27r and n is an integer number, which should not 
depend on the space coordinate ~r . The changes of A(/9 (7*) by discontinuous jumps means 
a divergence of its gradient, and as a consequence to zero contribution to the functional 
integral. Thus, the dependence on the constant integer number n can be included into the 
renormalization of the chemical potential /z, and the term 5ip/T~2'K /T can be neglected for 
large T. 

Integration over r/o in the generating functional < Z [J, V] >v finally gives 
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where 



< Z [J, V\ >y= \ DeDf exp {-%/ + iSj} 



(33) 



Seff = ^ [ d'r [L iV) + L (r)]' (34) 



and Sj is also expressed in terms of the fields e, /. The effective action Seff can be written 
in the form 

2~d/2 



%/ = ^(^fJ jd'^RLiffiR) (35) 
where Eq is the Larkin energy, which characterizes the scale of the fiuctuating random 
potential V . 



and Leff{R) is dimensionless representation of the function ^ (L (T*) + L (^)): 



L,jf(R) = ±1 - ( V^e) ( V^e) + / ( V«) ' / - ^nf (/)^ (36) 



2 

l/^l 

The sign of ±1 here is the sign of the chemical potential n, and R is the dimensionless 
coordinate variable 



can be fixed to unity if we choose Uc as 

,,(0) _ \A 

and the Lagrangian Lg// takes the form 



R = T'^MiA (37) 

(0) 



and the gradient operator is Vr = The dimensionless interaction constant a = -rin. 



n(°) = ^ (38) 
9 



L,ff(R) = ±1 - ( V^e) ( V^e) + / (v«)' / - {ff (39) 

Note that the effective action Eq. (153]) is derived without any assumptions about weakness 
of disorder and can be used for the description of the transition from strong localization to 
delocalization and superfluidity. 



IV. CALCULATION OF CHEMICAL POTENTIAL. 



In this section, as an application of the derived effective potential, we calculate the density 
of particles in three spacial dimensions {d = 3). The correlator < ^ (T*, t) \E' ("r", t) > 
defining the density of the particles, can be written in the form of functional integral 

<^(r,t)M/(r,t) >= j DeDeDfDj{f{-F')fexp{-Seff} (40) 

We consider the case of the large negative chemical potential /i < 0, | /i |>> £"0- The total 
number of particles is defined by this correlator 



N = y rfV < ^ (7", t) ^ ( r", t) > (41) 

and below we calculate the chemical potential as a function of the average particles density 
n = N/Q, where Q is the total space volume of the system. 

It is clear that only e, / field configurations with non diverging L^jf will contribute to the 
functional integral fj40l) . In particular, it means that Lefj (^-^^ should tend to zero rapidly 
enough for R 00. Thus, for large R, we have the equation 



- ( V^e) ( V«e) + / (Vr) / - 1 - iff = 

Solution of this equation defines the asymptotic R 00 of an instanton. In the whole 
region of distances R the instanton is defined by extremum of the action Seff- Lets find 
solutions with e = const. Then the equation takes the form 

7(v^,)V-i-(/)' = o 

Using the constraint // = 1, we obtain 

V«y/-/-(/)' = (42) 

We will look for solution in the class of spherically symmetric functions of the form / = 
exp(|A). Then Laplacian is 

Thus, for large R ^ 00 Eq. (H2l) can be written as 
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/ - / - ifr = (44) 



Moreover, we will look for solutions with the asymptotic 

/ (R) = exp (-i?) as R ^ oo (45) 
which means that the third term in Eq. fH4|) can be neglected at i? — > oo. The center of the 



instanton is an arbitrary parameter T^o? or its dimensionless analog Rq = A/2|/i|T^o- 



f i^Rj =exp\^-\R - Ro\j as R oo (46) 

and the integration over Rq is the integration over the zero mode. All instantons which 
differ only by parameter Rq give the same contribution to the functional integral, due to 
the translational invariance of the effective action. 
So we get in quasi classical limit for the correlator 



< ^(7',t)^(7',t) > (47) 
= J d^R, {r- (r) ) %xp I - (M) J d^R (L--) ' 

where 

j^inst -yinst ^'^^^ j-inst ^jinst^'^ 

One can look for the instanton configuration in the form 

exp (-\R - Ro\) 
finst (r\= V ^ (48) 



\R - i?oP + x2 

where the parameter x should be found from the minimization of the action Seff ( I35ll 
calculated on this configuration. We found that x = (f)^''^- The effective action Sejf for 
this value of x is 

1/2 



where 
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4 



(49) 



Finally we get the following estimation for the correlator 




2 



exp 




(50) 



Hence the total number of particles 



N = n(°) j dh j DRo (^r'' p)) exp 
= finlo)exp|-rMV^'Al 




) 



1/2 



A 



(51) 




Thus using Eqs. (15T|) and (l38ll we get the equation for the chemical potential, which is 



Strictly speaking, this equation is valid for the large values of the chemical potential so far 
as the saddle point calculation of the functional integral (HOl) is valid for >> E^. However, 
we think that, as one can guess from the form of the effective action (!35l) . the equation of 
this form is valid for |yu| ng Eq too. In this case the effective action is dimensionless 
and does not depend on any parameter. The main contribution to the functional integral 
f HU]) in this case provided by the field / with the asymptotic fHB]) . as as it is for instantons 
in the saddle point approximation. This contribution should be of the order of unity, so 
the estimation for the chemical potential is \fi\ ^ ng ^ Eq. Eq. (1521) has a solution for the 
densities smaller than the critical density n < n^., and does not have a solution for n > 
where nc ^ Formally, solution of Eq. (152|) . if exists, has two branches (see Fig.l for 
illustration). But the situation is different from the first order phase transitions, since only 
one of them is stable. The second solution has negative compressibility |^ < 0, and should 
be rejected. In a sense, the situation recalls a bifurcation point rather a first order phase 
transition. The stable solution of (|52|) with respect to as a function of n, in the region 

» Eq looks like 




(52) 



\ngj 
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(53) 




Fig.l Chemical potential function of density n. 



m, |27| 



It coincides with the expression for the chemical potential obtained recently in 
The value of the chemical potential for densities n < ric can be estimated as /i ~ fic, 
where /Xc ~ —Eq. When the density n becomes more than ric there is no solution of Eq. 
(132]) . which means that for n > the chemical potential will change the sign to positive 
and the system becomes delocalized and superfluid. 

Note that the exponential in Eq. fl5T]) is typical for the tail of the density of states in 3D 
for the negative energy of particles in a random potential. 

In order to understand long distance correlation properties of the system in the regime 
n < He consider the correlator < ^ {~ri^t) ^! {1^2, t) > for | T^i — 1^2 \^ oo. The functional 
integral in the form of the new variables is 

<^(T'i)^(T'2) >= j Df Deexp{-Seff} f (T'l) f {T'2)eiT' 1)6(1^ 2) 

In the saddle point approximation the configurations of fields (/ (^1) , e (V'l)), 
(/ ("^2) , e ("^^2)), when | "r'l — ^2 |— ^ 00, are described by the instantons localized near far 
distant points "r^i and 1^2, respectively, and the correlator < \E' (^1) ^ (^2) > is expone- 
tially small, i.e. the system has no long range order for densities n < ric. Note that the 
correlator < e (T^i) e (^^2) > is equal to zero for | T^i — T*2 |— ^ 00 due to the independence 
of the phases of the different instantons. 
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V. DISCUSSION 



The effective action for many-particle system in random potential is derived in the as- 
sumption of weak interaction between particles but without restrictions on the strength 
of the random field fiuctuations. The approach we develop is different from the usual Q- 
Lagrangian approach used in the theory of weak localization. The order parameter which 
could characterize the transition from localized to delocalized state has another structure 
and characterized by the "classical" field We obtain the effective action by the integration 
over the random external field at the first stage, following by integration over the quantum 
fiuctuations described by the field ip. The approximation we use for the integration over 
the field tjj is correct in quasi classical well as for small density of the system or 

weak interaction between particles. In this paper we consider the properties of Bose gas in 
random external potential with the Gaussian correlator for the random field fiuctuations in 
the case of the negative chemical potential, i.e., in the region of the strong localization. The 
obtained effective action has instantons which are localized in space and give the finite value 
for the effective action. The contribution of these instantons to the functional integral fHU]) 
defines the equation for the chemical potential as a function of the averaged density. The 
solution of this equation exists only if the averaged density is smaller than the critical value 
nc ~ The chemical potential for the densities n ~ ric is of the order ~ /^c = —Eq. 

We have shown that the system has no long range order for the densities n < Uc, whilst the 
chemical potential becoms positive, and the strong localized state transfers to the delocalized 
superfiuid state. 

One of the main important problems remaining for future is a more careful semiclassical 
investigation of instantonic solutions found here for different space dimensions. It would 
be also interesting to apply the developed method to the regime fi ~ /ic, where it could be 
possible to see the phase transition, probably of the first order, when fi changes the sign. 

Another problem where the described path integral Keldysh-Schwinger technique can be 
applied, is fermionic particles in a random potential. Here higher order powers of electron- 
electron interaction cannot be neglected, and the full action should be treated non per- 
turbatively. It seams reasonable to concentrate on the 1-1-1 dimensional space case for this 
problem and try to analyze a relation of the found chiral Lagrangian to solvable sigma model 
cases. 
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